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Abstract 

This paper shows how the enclosure method which was originally introduced 
for elliptic equations can be applied to inverse initial boundary value problems for 
parabolic equations. For the purpose a prototype of inverse initial boundary value 
problems whose governing equation is the heat equation is considered. An explicit 
method to extract an approximation of the value of the support function at a given 
direction of unknown discontinuity embedded in a heat conductive body from the 
temperature for a suitable heat flux on the lateral boundary for a fixed observation 
time is given. 
AMS: 35R30, 80A23 

KEY WORDS: inverse initial boundary value problem, heat equation, parabolic 
equation, thermal imaging, cavity, corrosion, enclosure method 

1 Introduction 

The aim of this paper is to show how the enclosure method which was originally introduced 
for elliptic equations in [3] can be applied to inverse initial boundary value problems for 
parabolic equations in multi-dimensions. We present only a simple case just for the 
explanation of the idea, however, the same idea will work also for more general cases. 

Let Q be a bounded domain of R m , m = 2, 3 with a smooth boundary. Let D be an 
open subset of Q with a smooth boundary and satisfy that: D C Q \ D is connected. 
We denote the unit outward normal vectors to dQ and dD by the same symbol v. Let T 
be an arbitrary fixed positive number. 

Given / = f(x,t), (x,t) G dQx ]0, T[ let u = u(x,t) be the solution of the initial 
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boundary value problem for the heat equation: 

d t u - Am = in (Q \ D) x ]0, T[, 

du 



dv 

du 
du 



on&Dx]0, T[, 
= f ondnx}0, T[, 



(1.1) 



u(x,0) =0 inQ\D. 
This paper considers the following problem for an explanation of the idea. 

Inverse Problem. Assume that D is unknown. Extract information about the location 
and shape of D from the temperature u{x, t) and heat flux f(x, t) for (x, t) G dVl x ]0, T[. 

The set D is a model of the union of unknown cavities or subregions where the heat 
conductivity is very low compared with that of the surrounded region Q\ D. Thus the 
problem is a mathematical formulation of a typical inverse problem in thermal imaging. 
The method to solve this inverse problem may have possibility of application in nonde- 
structive evaluation. 

There are extensive studies on the uniqueness and stability issues on this kind of 
problems. See [9] and references therein for the issues. In this paper we are concerned 
with the reconstruction issue. Note that in [5] Ikehata studied this type of problems 
in one-space dimensional case by using the enclosure method. In [6] he raised several 
questions related to the method in mult i- dimensional cases. However, it is still not clear 
that the method works also for two or three-space dimensional cases since the study in 
[5] fully makes use of the speciality of one-space dimension. 

Our main result is the following formula. 

Theorem 1.1. Given uj G S m ~ l let f be the function of (x,t) G dQx]0, T\ having a 
parameter r > defined by the equation 



dv 

f(x,t;r) = —(x)(p(t), 



(1.2) 



where v(x) = e^ x ' u ' and (p e L 2 (0, T) satisfying the condition: there exists /i G R such 
that 



lim inf r M 



f 

Jo 



-.-Tt 



<p(t)dt 



> 0. 



1.3) 



Let Uf — Uf(x,t) be the weak solution of (1.1) for f — f(x, t; r) and hoiyj) = sup^^, x-uj. 
Then the formula 



lim — — 



log 



an Jo 



-Tt 



dv \ 

v(x)f(x, t] r) — Uf(x, t) — (x) I dtdS 



h D (u), (1.4) 



is valid. 



Note that: if ip(t) is smooth on [0, T — 5[ with < 5 < T and t = is not a zero point 
with infinite order of (p(t), then (1.3) is satisfied for an appropriate /i > 0. 
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In [5] Ikehata considered the corresponding problem in one-space dimension for the 
equation d t u = iju x ) x with a smooth 7. Theorem 1.1 can be considered as an extension 
of (2) of Theorem 2.5 in [5] in the case 7 = 1 to two and three-space dimensional cases. 
Note that therein the parameter r 2 plays a role of r in Theorem 1.1. 

A brief outline of this paper is as follows. Theorem 1.1 is proved in Subsection 2.2. 
The proof is based on an integral identity which is described in Subsubsection 2.2.1. Using 
the identity, we give an asymptotic representation formula of the integral 

v(x)f(x,t;r) -u f (x,t)^(x) \ dtdS 

together with an estimate of the leading term which is an integral involving two Neumann- 
to-Dirichlet maps for the operator A — r. The estimate of the remainder term is described 
in Subsubsection 2.2.2 as Lemma 2.1 and that of the leading term is found in the same 
subsubsection as Lemma 2.2. To establish Lemma 2.1 we require an energy estimate 
of Uf(-,T): this together with the meaning of the weak solution of (1.1) is found in 
Subsection 2.1. In Section 3 we show how Theorem 1.1 can be extended also to extract 
another information about D, that is the distance from a given p outside Q to D. 



2 The enclosure method 

In this section firstly we specify what we mean by the weak solution of (1.1) and describe 
a necessary estimate. Secondly by giving the proof of Theorem 1.1 we show how the 
enclosure method can be applied to the heat equation. 



2.1 Preliminaries about the direct problem. 

We follow [1]. For two separable Hilbert spaces V and H with V C H C V, and 
a positive number T, the space W(0,T;V,V) is defined by W(0, T; V, V) = {u\u e 
L 2 (0,T; V),v! E L 2 (0,T; V')}. Note that v! means the derivative in t e]0, T[. 

Given / e L 2 (0, T; H~ l l 2 (d^l)) we say that u E W(0, T; H\n \ D), (H\n \ D))') 
satisfy 

d t u- Aw = in(fi\S)x]0, T[, 
du 

— = Oon0Dx]O,T[, (2.!) 

— = fondnx]0, T[ 

in the weak sense if the u satisfies 

< u'(t), > + / _ Vu(i, t) ■ V^(x)dx =< fit), V \ dn > in (0, T), (2.2) 
Jn\D 

in the sense of distribution on (0, T) for all ip G H 1 (il\D). 

Note that by Theorem 1 on p.473 in [1] we see that every u e W(0, T]H l (VL \ 
D), (if 1 (Q\D)) / ) is almost everywhere equal to a continuous function of [0, T] in L 2 (£l\D). 
Further, we have: 

W(0, T; H\Q \ D), (H\Q \ 33))') ^ C°([0, T];L 2 (n \ D)), (2.3) 
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the space C°([0, T];L 2 (Q \ D)) being equipped with the norm of uniform convergence. 
Thus one can consider u(0) and u(T) as elements of L 2 (Vt \ D). Then by Theorems 1 and 
2 on p. 512 and 513 in [1] we see that given u G L 2 (Q\ D) there exists a unique u such 
that u satisfies (2.1) in the weak sense and satisfies the initial condition u(0) = Uq. 

Let Mo — 0. By Remark 2 on p. 512 and Theorem 3 on p. 520 in [1] we have the 
continuity of w on /: there exists aCr>0 independent of / such that 

\\ u \\L2(o,T;m(n\D)) ^ CT||/||L2( 0) T;H-i/2(an))- (2.4) 
Moreover, from (2.2) and (2.4) we have 

IKIlL2( ,T;i?i(nVD)') - C7 1 ||/||L2( 0iT;// -i/2( 9n )). 

This together with (2.3) and (2.4) yields one of the important estimates in the enclosure 
method: 

\ u ( T )\l*(SI\D) ^ C T \\f\\ L 2^ jT . H ^i/2 i9n) y (2.5) 

In the following subsection we denote by Uf the weak solution of (2.1) with u(0) = and 
this is the meaning of the weak solution of (1.1). 

2.2 Proof of Theorem 1.1 

Define 

rT 

Wf(x;r) — / e~ Tt Uf(x,t)dt, x E Q 

J 

and 

g f (x;r) = / e~ Tt f(x,t)dt, x e dQ, 
Jo 

where r > is a parameter. This type of transform has been used in the study [5] for the 
corresponding problem in a one-space dimensional case. 

2.2.1 The basic identity 

The function Wf — w satisfies 

(A - t)w = e- TT u f (x, T) in Q. \ D, 

dw ^ 
— = ondD, 
ou 

— — = Of on oil. 

du r 

Let v — v(x) satisfy (A — r)v = in Q. Integration by parts yields 

f ( Q fV — w f^—\ dS = — [ Wf—dS + e~ TT [ Uf(x,T)v(x)dx. (2.6) 
Jan \ du J J 3D du Jn\D 
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Let pf = p be the unique solution of the boundary value problem: 

(A-r)p = inQ\D, 

^ = ondD, 
dv 

— — = Of ond\l. 

dv y/ 

Set €f — Wf — pf. Since we have 

L { 9fv ~ pf t) dS = - f B D pf t dS > 

from (2.6) we obtain the basic identity: 

L { 9fv ~ wf t) dS = L { 9fv - pf t) dS 

dv 

d£> + e ' " 

ln\D 

Note that e/ = e satisfies 

(A - r)e = e- rT M/(x, T) in Q \ D, 



f dv f 
— / €f—- dS + e~ rT / Uf(x,T)v(x)dx. 
J 3D dv Jn\D 



(2.7) 



de 

dv~ = 0ondD > (2-8) 
de 

tt = on 0ft. 



2.2.2 Two lemmas 

Now we choose a special / having the form: 

dv 

f( x ,t;r) = —(xMt), (x,t)edQx]0T[, (2.9) 

where <p G L 2 (0, T). 
Since 

dv fT 
9fi x \ T ) — ~q^( x ) J e- Tt ^{t)dt, 

we have the expression of the first term of the right hand side of (2.7): 

l fa - «s is = r /„ i - ^t ds - i2M) 

where Rq(t) and Rd{t) are the Neumann-to-Dirichlet maps on dVl for the operator A — r 
in ft and Vt\D with the homogeneous Neumann boundary condition on dD, respectively. 
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Hereafter we choose a special v in (2.9) which is a solution of the equation (A — t)v = 
in R m : 

v ( x ) = e^", 

where uj G S" 1 ^ 1 and r > 0. Thus, we choose f(x,t; r) as described in (1.2). 

The following is an easy consequence of the estimate (2.5), an a priori estimate for 
(2.8) and the growth order of v as r — > oo. 

Lemma 2.1. We have, as r — > oo 

- I e f ^dS + e' rT [ u f (x, T)v(x)dx = Oi^ e~ rT e 2 ^ h ^) 

J 3D OV Jfl\D 

where 7 is a constant independent of r and h^{uj) = sup^^rr • uj. 

The integral in (1.4) is just the left-hand side of (2.7). This fact, (2.10) and Lemma 
2.1 imply that 

j m £e- Tt (v(x)f(x,t;r) - Uf (x,t)^-(x)\ dtdS 

(2.11) 

= [ T e-^(t)dt[ ^(R,(r)-R D (T))^dS + 0(r' e ^ T e 2 ^ h ^y). 
Jo Jan ov ov 

Thus, the main contribution on the limit (1.4) is the term given by two Neumann-to- 
Dirichlet maps. For this term we need the following estimates from the both sides: 

Lemma 2.2. There exist /i l5 ji 2 G R, r > 0, Ci > and C 2 > such that, for all r > r 



dv 
an dv 



R D (r))^-dS 



< C 2 T>». 



From Lemma 2.2, (2.11) and (1.3), we have formula (1.4) in Theorem 1.1. Note that 
in Lemma 2.2 the lower estimate is essential and the strict values of n\ and /i 2 are not 
important to obtain formula (1.4). 

The proof of Lemma 2.2 given below follows the argument in the enclosure method 
applied to the Hemlmholtz equation (A + k 2 )u = [3, 4]. It became much simpler 
compared with that for the Helmholtz equation since our operator A — r is strictly 
negative for r > 0. 

Proof of Lemma 2.2. Consider the solution R(x) of the following boundary value problem: 

(A-r)i? = inQ\D, 

OR dv An 
ov ov 

^=0ondn. 
Ov 

Then we have, given r\ > 

\ R \m(n\D) < C(\Vv\ L 2 {D) + t\v\ L 2 (d) ) 
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where C > is independent of r > rj. Noting that R4(r)(dv / 'dv) = v, RD{r)(dv / dv) = 
v — R on dfl and the functions R and v are real-valued, we have the integral identity 

f dv dv 

= I _\VR\ 2 dx + r ( _\R\ 2 dx+ [ \Vv\ 2 dx + r [ \v\ 2 dx. 
Jn\D Jn\D Jd Jd 

A combination of them gives 

2r / e 2 ^ x ^dx < [ ^(R$(r) - R D (r))^dS < C 2 t 3 [ e^^dx. 
Jd Jan dv dv Jd 

Now the conclusion follows from the corresponding estimate for the integral of e 2 ^'^ 

over D (cf. Propositions 3.1 and 3.2 in [3]). 

□ 



3 The use of another v 

It is possible to use different v from the exponential solution. In [7] it is shown that: 
given p G R m \ Q one can construct a«G H 2 (fl) depending on r > and satisfying the 
equation (A — r)v = in Q such that 

v (x) = e -^- p l j|x - p\-^-Vn + o | . 

Note that the leading term never vanish on Q; in the case when m = 3, one can drop the 
term 0{l/^/r) in the v above. 

Using this v, we obtain the following formula. 

Theorem 3.1. Let p G R m \ and replace v of f in (1-2) with the v above. Let 
Uf = Uf(x,t) be the weak solution of (1.1) for this f = f(x,t;r). Then assuming (1-3), 
one has the formula 



lim — — log 

— >™ 2^ B 



[ [ T e~ Tt (v(x)f(x,t;r) - u f (x, t)^-(x) ) dtdS 
Jan Jo \ dv 



= ~d D (p), 



where dr>(p) denotes the distance from p to D, 

d D (p) = inf{|y - p\ \ y E D}. 

Proceeding once more by the method used to prove (1.4), one knows that the key of 
the proof of Theorem 3.1 is the following lower estimate of the integral of e~ 2 ^ x ~ p \ over 
D. 

Proposition 3.2. There exists fi G R such that 

liminf rV^ dc(p) / e'^^dx > 0. (3.1) 

T >00 J D 
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Proof. Choose x G 3D such that do(p) — \ x o — p\- Since dD is smooth, one can find 
an open ball B with x G dB such that B C D. Since the integrand is nonnegative, it 
suffices to prove (3.1) in the case when D = B. We observe that 

| Ve -2 > /T|x-p|| = 2^fre- 2 ^ x - p y 
Using the co-area formula (Theorem 2.7.1. on page 76 in [10]), we obtain 




x G B | e - 2 ^-Pl = t})dt 



(3.2) 

POO 

= / H m ~\{x eB\\x-p\ = s})e- 2 ^ s ds, 

Jd D {p) 

where H m ~ l denotes the m — 1 dimensional Housdorff measure. Since for appropriate 
H' > and C > H m ~\{x £ B \ \x — p\ = s}) > C(s - ^(p))* 4 ' as s j d D (p), from (3.2) 
we obtain the desired estimate. 
□ 

4 Conclusion and Remarks 

The procedure of extracting the support function of D is extremely simple and summarized 
as follows. 

(i) Give the direction oo G S' m_1 . Fix a large r > and give the heat flux across dVt over 
the time interval ]0, T[: 

dv 

f(x,t;r) = fo{x)<p(t), (x,t) G dQx}0, T[, 
where v(x) = e^ x ' UJ and ip(t) satisfies (1.3) for a //. 

(ii) Measure the temperature Uf(x,t) on dfl over the time interval ]0, T[. 

(iii) Compute the quantity 

/ [ T e~ Tt (v(x)f(x,t;r) - u f (x, t)^-(x) ) c?to?5 

JO \ OV J 

as an approximation of h D {uj). 

As a corollary of Theorem 1.1 we have a constructive proof of the uniqueness of re- 
covering the convex hull of D from infinitely many sets of the temperature and heat flux 
on dVt over the time interval ]0, T[. Thus our result can be considered as an extension of 
the enclosure method of infinitely many measurements version [3]. 

The method will cover also more general cases without serious difficulty: inclusion, 
parabolic equations with variable coefficients, Robin boundary condition. And it may be 
possible to apply our method to the corresponding problems for hyperbolic equations and 
systems, Stokes system, some governing equations in thermoelasticity, Maxwell systems, 
etc. Such cases will be reported in detail in forthcoming papers. 

Note that in [8] we have already applied the enclosure method of a single measurement 
version [2] to the problem in three-dimensions. It means that therein the heat flux / is 




8 



fixed and independent of large parameter r. The information extracted therein is different 
from the information obtained in this paper. The analysis developed therein is based on 
the potential theory and quite delicate. See also [7] for an application of the enclosure 
method to an inverse source problem for the heat equation in multi-space dimensions. 
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